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ABSTRACT 

Based on the earlier work of Gunn and McCrea we model the formation of 
globular clusters in merging galaxies. Neutral hydrogen observations of dwarf 
irregular galaxies as well as more luminous systems are used to provide the key 
parameters of the model. The observations indicate that clusters with the mass 
of globular clusters should still be forming today. The model is incorporated into 
a phenomenological picture of galaxy evolution making use of a simple chemi- 
cal evolution model. These results are compared to recent observations of the 
metallicity distributions of F and G stars from a recent large SDSS survey. The 
comparisons are consistent with an anisotropic collapse and merging of a large 
number of dwarf irregular galaxies for the formation of the Galaxy. 

Subject headings: Galaxy: evolution — Galaxy: halo — globular clusters: general 



1. Introduction 

One of the remarkable observational facts about the stellar content of galaxies is the 
ubiquity of globular star clusters. These clusters can be found in almost all but the least 
massive galaxies. Generally globular clusters come in two major families of comparable 
mass: metal poor, old blue clusters and metal richer, younger, red clusters. A comprehensive 
review of globular cluster systems is given by Harris (2001). Below we present quantitative 
arguments for the formation of clusters during mergers and use the observations to obtain the 
key parameters of the model. This formation picture leads naturally to a phenomenological 
model for galaxy evolution. Combined with a simple chemical evolution model comparisons 
with recent data for Galactic halo stars are made. 
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2. Merger Induced Globular Cluster Formation 

The possibility of forming globular clusters from collisions between gas rich bodies was 
considered by Gunn (1980) and independently by McCrea (1982). We revive and extend the 
Gunn and McCrea analysis here. 



2.1. The model 

Within each of two equal mass interacting galaxies consider an undisturbed column of 
gas of length 1, cross-sectional area A and density po meeting its counter- part with relative 
velocity 2V (see Fig.l). During the collision, the gas is compressed to column length A and 
density p. The shocked gas is assumed to cool to T~ 10'' degrees. This assumption will 
be justified later-suffice it to say that the cooling time within the compressed volume {tcooi) 
must be shorter than the collision time {tcou ~ 1/^) for the assumption to be valid. By 
conservation of mass and for simplicity assuming A to remain constant 

2polA = pXA (1) 

equivalently in terms of the column density 

2E„ = E (2) 

and 

Po A 

with the relative velocity being 2V, the relation between po and p becomes (e.g. Spitzer 
(1978), equation 10-24) 

pAV' ^ p (^) (4) 

so that 

— ~ 0.048 V^/(TJii) (5) 

po 

where the bracketed quantity in equation (4) represents the square of the sound speed, and 
in equation (5) the unit of V is km/sec, T4 is the temperature in units of 10^ degrees and /j 
is the mean molecular weight of the shocked gas. The quantity (T4//i) occurs frequently in 
the expressions below. With p, — 0.59 a typical value for this quantity is ~ 1.7. Combining 
(3) and (5) 

J ~ 0.024 Vy{T,/p) (6) 
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The Jean's length, Aj is 



/ nkT V'^ 246 



with p in units of M^jpc" . Equating A with Aj eliminating p 

6.05 X lOS^ , , 
A = Aj = ^^(r4^)pc (8) 

Finally, the Jean's mass is given by 

Mj = p}? = 2EoA2 = 6.05 X 10^A(r4///) Mq (9) 



If the cooling time of this gravitationally unstable mass, Tcooi-, is shorter than the free-fall 
time (t//) then fragmentation can occur. The cooling function A(Z) which when multiplied 
by the square of the particle density of hydrogen gives the cooling rate per unit volume, falls 
precipitously from a metallicity independent peak near T~T4 to a metallicity dependent 
plateau below (e.g. Dalgarno & McCray 1972). The coohng rate from heavy elements alone 
at r4 ~ 0.9 is > 4 X 10-26(^/^^)^2^ (e.g. Fall & Rees 1985, §IV, p. 22). From t^ooI = 
1 .5nkT / (An'jj) and with the density in M© pc~^ the cooling time scale becomes ~ 6.96 x 
10^{Z/Zq)-\T^/ii)/p yrs. Similarly t// = (32Gp/37r)-V2 reduces to ~ 8.08 x IOVp^^^ yrs. 
Finally, from (3) and (4) one can write Tcou = t^S? where = X/cg ~ 1.07 x 10^\{T4/ p)~^^^ 
yrs, Cg is the sound speed, A is in pc, and 3? is the Mach number {2V/cs) and is ^ 1. 
Eliminating the density using (7) the relevant time scale ratios become 

Tcool/Tcoll < Tcooi/Ts ~ 1-07 X lO'' X(Z / ZQy\n/ p)'/'' (10) 

and 

Tcooi/rff ~ 3.5 X WX{Z/ZQ)-\T,/py/'' (11) 

As will be shown in the next section, the observationally determined value of A is such that 
both of these ratios are less than unity. Hence our assumption that the shock is isothermal 
is justified. We are also assuming that the temperature of the gas remains constant at 
the maximum of the cooling function near T4 ~ 1 until the mass becomes Jeans unstable. 
Further, one can expect fragmentation to occur when the gravitationally unstable cloud 
begins to collapse as long as Z/Zq is not much less than ~ 0.01. This is also the threshold 
metallicity below which globular clusters are not observed. 

Assuming that stars are eventually formed, another consideration is the efficiency 77 with 
which gas is turned into stars. This is not well known at present so for the later discussion we 
assume this quantity to be 77 ~ ^star/^baryon ~ 0.07. This ratio represents the global stellar 
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mass density divided by the total baryon density. A further question then arises because 
the mass loss associated with star formation and evolution can lead to the disolution of any 
cluster of stars formed. This problem has been investigated by Geyer & Burkert (2001) 
and more recently in their discussion of the globular cluster mass function by Parmentier & 
Gilmore (2007). Unfortunately the observations which would allow an empirical estimate of 
the bound to unbound stellar ratio are not yet available so that more sophisticated modelling 
is required to answer this question in the context here. To proceed we shall consider rj to set 
an upper limit to the mass of a bound cluster of stars formed. 

A final point to note is that in this model the gas is concentrated by collision and the 
co-orbiting stars and dark matter will at least initially be unaffected by this process. Hence 
we do not expect the resulting globular clusters to contain a significant amount of dark 
matter. 



2.2. Relating the model to the Observations 

It is clear from the above discussion that the crucial quantity to determine is A while the 
principal variables are a characteristic length 1, and a characteristic velocity V. Given that 
the available observational parameters of the dominant gaseous component of galaxies are 
the radius of the HI disk, ichi, the mass of the HI gas, M^/, and the peak rotational velocity 
of the HI, Yrot, we begin by identifying 1 with thi and V with Yrot- By identifying V with Yrot 
wc arc considering mainly edge-on collisions. While the mass does not directly enter into the 
discussion above, correlations between it and the other variables lead to interesting results. 
The mass parameter actually used here is I.SMhi in order to allow for the contribution of 
helium. Throughout the rest of this work, the symbol Mfji will generally represent both the 
hydrogen and helium gas content. 

The data discussed here comes from two sources. The first for relatively bright galaxies 
are from Broeils and Rhee (1997) and the second for dwarf irregular galaxies are from Begum 
et al.(2008a,b). In order to intercompare both data sets a correction to both r and M was 
made for the small difference in assumed Hubble constant. In addition because the HI radii 
were measured at a different limiting surface brightness, a correction of -1-0.178 in log r was 
applied to the Broeils and Rhee data. This was determined by jointly correlating the radii 
and the masses and finding the minimum dispersion. This correction is well within the 
uncertainties quoted by Begum et al. in their discussion of this problem. Fig. 2 shows the 
resulting relation between 2 log yhi and log Mhi- The solid line has a slope of 1 while a least 
squares fit for this slope gives 0.994 ± 0.012 for all 137 galaxies. The logarithmic intercept 
(log M-21og r) is 0.846 ± 0.132. Hence in agreement with Begum et al. the HI surface 
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density appears to be remarkably constant over nearly 5 orders of magnitude in mass. Note 
that this face-on surface density is not to be confused with E,, above where it represents the 
column density. 

Equation (6) above implies that if A is a constant then a relation between 2 log V and 
log r should determine its value. Fig. 3 shows such a plot. Possibly a more famihar way to 
look at this relation is to eliminate r by using the relationship established in Fig. 2. Since 
roc M^/^ equation (6) implies Moc V^^^ i.e. an HI Tully-Fisher relation is predicted. This 
is shown in Fig. 4. The solid line in both Figs 3 & 4 has a slope of 1. Least squares fit 
to the data with Vrot > 41.7 km s^^ (i.e. above the dashed lines) give 0.978 ± 0.078 in 
Fig. 3 and 0.989 ± 0.089 in Fig. 4. Both regressions included 109 galaxies. The logarithmic 
intercepts are (log r-2 log Vrot)^ 0.168 ± 0.250 and (log M-4 log V)=1.18 ± 0.515. Solving 
for A from each one obtains A = 61.3 ± 35.3(T4///) pc and A = 61.4 ± 36.3(T4//x) pc 
respectively. Although the uncertainties are large the numerical difference is small given the 
tight correlation seen in Fig. 2. Given this value of A the time scale arguments above are 
justified. It should be noted however that the condition for fragmentation will be violated for 
values of Z/Z© not much below 0.01. This provides a natural explanation for the observed 
abundance threshold of this order below which globular clusters are not observed a point 
also made by Gunn (1980). 

Returning to equation (9) the Jean's mass becomes Mj = 3.7 ± 2.1 x 10^(T4/^)^ Mq. 
Allowing for a gas to star formation efficiency of r/ ~ 0.07 results in a predicted globular 
cluster mass upper limit of Mqc ~ 2.6 ± 1.5 x 10^(r4//x)^ Mq. 

Begum et al. (2008a) have also noted that the dwarf irregular galaxies with the lowest 
rotational velocities (i.e. those below the dashed line in Figs 3 & 4) lie below the extrapolated 
baryonic Tully-Fisher relation. If all of the galaxies are rotationally supported, then V^^^ = 
aGMui/rHi where a is the ratio of the total mass enclosed within r^/ to the gas mass. Given 
the relatively tight relation in Fig. 2 the increased scatter in Figs 3 & 4 can be understood 
as due to variations in this ratio. In order to account for the data below the dashed lines, 
either the gas is not fully rotationally supported or there is a systematic shift in the value 
of a towards a lower value inspite of the tight correlation between the HI radius and mass 
in fig. 2. 

Suppose that the galaxies below the dashed lines are indeed relatively gas rich (i.e. 
low a). These objects could then be more representative of the universe as it was 13 
Gyr ago and be the proto-galactic clumps from which the present galaxies were assembled. 
Fitting lines of slope unity through these 28 points in figs 3 and 4 then yields values of A 
of 310 ± 212{T/^/ ji) pc and 307 ± 271(T4/yu) pc repectively i.e. a factor of ~5 higher 
than from the more massive galaxies and hence a Jean's mass higher by the same amount. 
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In order to recover the observed globular cluster masses one would need to invoke a star 
formation efficiency smaller by the same amount. In any event because the clusters are old 
and relatively metal poor, it will be argued below that these clusters were formed from the 
mergers of dwarf irregular galaxies whether those below the dashed lines in figs 3 and 4 or 
at the low mass end of those above. 

At the beginning of this discussion we chose to identify the 1 and V of the model with thi 
and Vrot of the gas rich galaxies. The choice of Vrot would seem to be natural since during 
the interaction between two counter rotating galaxies the undisturbed gas is rotationally 
supported. The choice of thi for 1 then leads to a 'reasonable' mass for a globular cluster. 
A further consistency check comes from the model prediction that the column density of 
gas (So = Pol) must be a constant and equal to = 3.02 x 10^(T4//i)/A. With the value 
of A derived above (assuming I—ichi) Eq = 496 ± 285 MqPc~^. This in turn imphes a 
value for the undisturbed particle density varying from ~ 7 — 0.2 gas particles/cm^ as the 
galaxy radius increases from 2.5 to 100 kpc. For a Milky Way sized galaxy, {thi 25kpc), 
n^ ~ 0.6 ± 0.4 cm~^ compared to an observed value of uh — 0.57 cm~^ (Dickey & Lockman 
1990). 

If structure in the universe forms hierarchically as is currently believed i.e. low mass 
structures form first and merge to make larger objects, then we expect the oldest clusters 
to have been formed at the time when the low mass systems were merging. In fact genuine 
globular clusters (both blue and red) arc found in many 'surviving' dwarf irregular galaxies 
(e.g. van den Bergh 2006, Gcorgicv et al. 2008). Generally, as larger structures merge, 
we expect clusters of similar mass to form but have progressively younger ages and higher 
metallicities. Today, given the right circumstances, these objects should still be forming and 
may resemble the super star clusters first observed by Schweizer (1987) in the Antennae 
galaxies. Conversely, the Galactic building blocks could not have been much more massive 
than the dwarf irregular galaxies before merging or the red globular cluster backbone would 
have been younger and more metal rich than observed. 

In the next section an outline of galaxy formation based on this picture is presented. 

3. Implications for Galctxy Evolution 

We now incorporate the above cluster formation model into an outline for galaxy evo- 
lution. The two separate globular cluster families (i.e. the old metal poor blue ones and the 
slightly younger metal rich red ones) are also found in different locations in a galaxy. Unlike 
the red clusters which are almost exclusively found in the inner regions of galaxy halos. 
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the blue globular clusters can be found in both the inner and outer regions (e.g. Harris, 
2001). In order to account for these differences we invoke an anisotropic collapse picture. 
We imagine that a spiral galaxy forms within a filament of proto-galactic clumps of gas and 
dark matter. The filament is assumed to collapse first in a direction perpendicular to the 
filament length and then along its length. Generally the angular momentum content of a 
galaxy is distributed so we expect the low angular momentum building blocks to collide in 
the first collapse and in a major burst of star formation form the blue clusters and associ- 
ated field stars. While the previously formed stars and dark matter will form a more diffuse 
background, the gas left over after the star formation burst from which the blue globulars 
and associated field stars form will fall to the center and eventually be available to form the 
nucleus of the galaxy and contribute to the low angular momentum component of the bulge. 
The end result is a population of field stars formed prior to the mergers and their associated 
dark matter halos and post-merger blue globulars and stars formed as a result of the merging 
activity but with different kinematics than the prior population. The higher angular mo- 
mentum building blocks continue to become enriched and slowly spiral in towards the center 
at which time they begin colliding with each other and leave behind a diffuse background of 
stars, dark matter and any previously formed blue globular clusters. The gas undergoes a 
second major burst of star formation and red globular cluster formation before the residual 
gas eventually falls to form the disk. A chemical cvohitionary model based on this picture 
which was also able to account for the global star formation history was given in H04. Here 
we compare the model to new and previously unavailable data based on a large survey of 
stars with photometric metallicity determinations from the Sloan Digital Sky Survey (Ivezic 
et al. 2008). 

Before discussing the observations we shall describe the chemical evolution model used 
to make the comparisons. It is similar to the model described in H04. There the equations 
are modified slightly from a one zone mass loss model in order to mimic a multi-zone model 
which works as follows. Suppose that in the one zone model the chemical evolution is 
interrupted before all of the gas has turned into stars as might occur during a merger for 
example. We are then left with a sharply truncated metallicity distribution of stars and a 
mass of gas with a single metallicity. Now suppose that these truncations occur over a finite 
range of metallicity dictated by a Gaussian distribution for example. Then we shall end up 
with a more smoothly truncated stellar distribution and finite spread in the metallicity of the 
left over gas which should be almost Gaussian shaped. The end result mimics the result of 
many zones having undergone truncations (mergers) over a range of metallicity. The model 
also allows for mass loss to occur at a rate proportional to the star formation rate by defining 
an effective yield. 

Fig. 5 is taken from the paper of Ivezic et al. (2008) (bottom right of their fig. 7). 
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The points with error bars show the metaUicity distribution of halo stars with height above 
the plane | Z |=5-7 kpc. The blue curve shows the predicted metallicity distribution of 
the left over gas after the collisions have taken place. Because the model does not predict 
either the metallicity distribution or the mass of the stars formed during the mergers we 
have assumed that the post-merger stellar and cluster distribution is the same as that of the 
left over gas and used a gas to star formation efficiency of ~ 0.07 to estimate the mass of the 
stars and clusters. The sohd black line (arbitrary normahzation) shows the distribution of 
metallicity of the stars formed prior to the mergers. It was calculated using exactly the same 
chemical evolution model given in H04. The model parameters are log(pe///^0) = —0.3, 
log(Zc/Z0) = —1.50, and Ouue — 0-32. Using the above star formation efficiency, the ratio of 
the mass of stars formed during the mergers to those formed prior to the collisions is ~ 0.8. 
This number will also depend on the value assumed for Pe//- For illustrative purposes here 
and to minimize parameters we have kept the effective yield the same as the true yield (i.e. 
O.5Z0). The stars formed prior to the merger are clearly more metal poor than the post- 
merger population and apparently do not show up in the data in Fig. 5 although no attempt 
was made to fit more than one component to the data. These stars are expected to have 
different kinematic properties since they would not have been affected in the same way by 
the collisions which produced the burst of star formation and clusters. In fact, these stars 
may belong to the second blue halo population which is more metal poor and with different 
kinematics identified in a deeper survey by CaroUo et al. (2007). 

Note that a prediction of this model is that generally the stars enclosed by the blue 
curve should have the same kinematics and spatial distribution as the blue (halo) globular 
clusters since they are assumed to have formed at the same time. Further these clusters 
should belong to the 'young' halo group of Zinn (1993). 

Fig. 6 from the top right of fig 7 of Ivezic et al. shows a similar plot but now with stars 
from I Z 1= 1.5 — 2.0 kpc. Here the red curve is the chemical model metallicity prediction 
assumed for stars and red clusters formed during the second set of collisions. Instead of using 
a Gaussian to determine the metallicity distribution of the stars/clusters formed during the 
collisions, we have used an extreme value distribution function which has the shape of an 
asymmetrical Gaussian, i.e. the function f of H04 is replaced by 

/ = exp{—exp{—{log Zc — log Z)/w)) (12) 

and the Gaussian function in equation (8) of H04 (df/dlogZ) is replaced by the expression 

df /dlog Z — exp{{—{log Z^. — log Z)/w) — exp{—{log Z^ — log Z)l/w))/w (13) 

The model parameters are log(pe///^0) = —0.3, \og{Zc/ZQ) — —0.68, and w — 0.16. This 
function could have been used for the blue clusters as well but the fit was marginally better 
using the Gaussian for them given that only one component was being fit. 
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The solid black line (arbitrary normalization) shows the metallicity distribution of all 
stars formed in the surviving population of building blocks prior to the second major merger 
phase. This group of stars (called the red spheroid in H04) will be more diffuse, have different 
kinematics and will possess a metal weak tail. The mass ratio between these two components 
(post to prior) for the parameters assumed above is ~ 0.2. 

Note that any blue clusters that formed prior to the first major merger phase and which 
remained bound to its lower mass galaxy while avoiding the first mergers will also have the 
kinematics of the red spheroid population and hence would be members of the 'old' halo 
component of Zinn (1993). The red clusters should have the same kinematics as the stars 
under the red curve of fig 6. The closest analog of a thick disk population in this model are 
these same stars. 

The main purpose of this discussion is to show that the phenomenological model outlined 
above is at least qualitatively consistent with the new data and with our picture of globular 
cluster formation. Further, the general outline for galaxy evolution described above may 
help explain additional globular cluster related questions. For example, the blue-red cluster 
phenomenom is also present in elliptical galaxies even though they are morphologically quite 
different from the gas rich systems. Could it be that eUipticals formed similarly but at the 
intersection of more than one filament? The increased frequency of collisions would convert 
more gas to stars and lead to a higher specific frequency of globular clusters. In both cases 
variations in the distribution of angular momcmtum content could explain the differences in 
the relative frequency of blue to red clusters. We leave these interesting questions to future 
work. 



4. Conclusions 

We have given a quantitative model for the formation of globular clusters based on the 
early work of Gunn & McCrea. Clusters are assumed to form in merger induced collsions. 
We make use of recent HI data for dwarf irregular galaxies and earlier observational work 
to derive the key parameters of the model. We then combined this model with a simple 
chemical evolutionary model to predict the metallicity distributions of stars formed before 
and during these mergers. These predictions were compared with relatively recent survey 
data and shown to be consistent with both models. As additional data becomes available it 
will be possible to refine these ideas into a more sophisticated picture for galaxy evolution. 

The author is grateful to Zeljko Ivezic for promptly sending the data contained in Fig. 
7 of their important paper. He also wishes to thank Brent TuUy for suggesting references to 
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the literature, Ray Carlberg for helpful comments and the referee for a constructive report. 
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Fig. 1. — Two gas columns approach one another with relative velocity 2V and coUide 
reducing the column length from 1 to A and increasing the density from po to p 
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log (1.3 Mf^,) 

Fig. 2. — Twice the log of the radial extent of the neutral hydrogen gas thi versus the log 
of 1.3x Mhi- The sohd line has a slope of unity. The sohd points are from Broeils & Rhee 
(1997) while the open circles are from the work of Begum et al. (2008a,b). The units of ichi 
and Mhi are pc and Mq respectively. 



5 




Fig. 3. — Twice the log of the rotation velocity of the neutral hydrogen versus the log of its 
radial extent. The sohd line has a slope of unity and the dashed hne is the adopted lower 
limit of the velocity below which the predicted relation is assumed to break down. The 
symbols are as defined in fig. 2. The unit of rotational velocity is km/sec 
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Fig. 4. — Four times the log of the rotation velocity versus the log of the mass of gas. The 
solid line has a slope of unity and the dashed line is the adopted lower hmit to the velocities 
of galaxies assumed to be obeying the HI TuUy-Fisher relation. The symbols are as defined 
in fig.2. 
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Fig. 5. — The data points are from fig 7 (bottom right) of Ivezic et al. (2008) and show the 
metaUicity distribution of stars within the distance interval above the plane of | Z |= 5 — 7 
kpc. The blue line is derived from the chemical evolution model under the assumption that 
the stars were formed in a major merger phase associated with the formation of the blue 
globular clusters. The black line shows the predicted distribution of metallicities of the stars 
which formed prior to the merger phase. The normalization of this curve is arbitrary. 
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Fig. 6. — The data points are from fig 7 (top right) of Ivezic et ah (2008) and show the 
metalhcity distribution of stars within the distance interval above the plane of | Z |= 1.5—2.0 
kpc. The red line is derived from the chemical evolution model under the assumption that 
the stars were formed in a second major merger phase associated with the formation of 
the red globular clusters. As above the black line with arbitrary normalization shows the 
predicted distribution of metallicities of stars formed prior to the mergers. 



